Theory of Underdoped Cuprates 



in 

ON 
On 

S3 
33 
in 

in 



> ■ 
in ■ 

o . 
\o • 
o ■ 
in : 

On . 
■ 

c3 ■ 



C3 
O 

o 



13 



Xiao-Gang Wen, and Patrick A. Lee 
Department of Physics, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139 

(June, 1995) 

We develop a slave-boson theory for the t-J model at finite doping which respects an 517(2) 
symmetry - a symmetry previously known to be important at half filling. The mean field phase 
diagram is found to be consistent with the phases observed in the cuprate superconductors, which 
contains d-wave superconductor, spin gap, strange metal, and Fermi liquid phases. The spin gap 
phase is best understood as the staggered flux phase, which is nevertheless translationally invariant 
for physical quantities. The electron spectral function shows small Fermi pockets at low doping 
which continuously evolve into the large Fermi surface at high doping concentrations. 

PACS numbers: 74.25.Jb,79.60.-i,71.27.+a 



The parent compound of the cuprate superconductors 
is an antiferromagnetic (AF) insulator. With hole dop- 
ing, AF is rapidly destroyed and a metallic state emerges. 
It is well established that at optimal doping, a Fermi sur- 
face exists with area 1 — x where x is the concentration 
of doped holes. [Q On the other hand, for x <C 1, the sys- 
tem remains AF with a doubled unit cell, and the Fermi 
surface is expected to be small pockets centered around 
the (§,§) point. || An important question is how the 
low-lying electron state evolves from small x to optimal 
doping. This intermediate region, called the underdoped 
regime, also exhibits unusual magnetic properties often 
referred to as the spin gap behavior. Unlike optimally 
doped systems, where the magnetic susceptibility \ is 
temperature independent, underdoped cuprates gener- 
ally show a reduction in \ at temperatures below 400K 
or so. H Below 150K, x an d the NMR relaxation rate 
decreases abruptly in an activated manner. It has been 
argued that this is observed only in bilayer material. ||] 
In this paper we shall concentrate only on the high tem- 
perature spin gap behavior, which we view as evidence 
for the formation of spin singlets within the Cu-O layer. 
We shall address the question of how the spin gap man- 
ifests itself in the electronic spectral function, and how 
the Fermi surface evolves from small pockets to a large 
Fermi surface which satisfies Luttinger's theorem. This 
last question was addressed in a weak coupling theory 
H where fluctuating spin density waves induce shadow 
bands. We would like to study the strong correlation 
limit which we believe to be more appropriate for the 
cuprates, and we take the t-J model as our starting point. 

A standard way of enforcing the constraint of no dou- 
ble occupancy of the t-J model is to write the electron 
operator c a i in terms of auxiliary fermions and boson 
particles c a i = f a ib\ and demanding that each site is oc- 
cupied by either a fermion or a boson. In a mean field 
(MF) treatment, the order parameters Xij — (faifaj) 
and Ay = (fufej — fnfij) describes the formation of 
singlets envisioned in Anderson's resonating valence bond 
(RVB) picture. || At zero doping, the translationally in- 



variant solution can be described as a 7r-fl.ux phase |J or 
a d-wave pairing state with |Ay| = \Xij\- The symmetry 
which underlies the degenerate MF states has been iden- 
tified as SU(2), which expresses the idea that a physical 
up spin can be viewed as either the presence of an up spin 
or the absence of a down spin fermion. § In the conven- 
tional theory, the SU(2) is broken to 17(1) upon doping, 
and only the d-wave state survives as the MF solution. 
||[)| This scenario has been used as an explanation of the 
spin gap phenomenon. || 

In this paper we present a new formulation of the con- 
straint which preserves SU{2) symmetry away from half- 
filling. Our hope is that since SU{2) is an exact sym- 
metry at half-filling, the MF approximation of the new 
formulation may capture more accurately the low energy 
degrees of freedom and may be a better starting point 
for small x. We are also motivated by the photo-emission 
experiment on the insulating cuprate, JTo[ ] which finds a 
large excitation energy at the (0,7r) point, comparable 
to that at (0, 0). This is just what is expected from the 
7r-flux phase spectrum, suggesting that the AF state may 
resemble the 7r-flux phase at short distances. As we shall 
see, the SU(2) formulation provides a scenario for how 
the 7r-flux phase is connected to the spin gap phase and 
how the hole pockets evolve upon doping. 

The SU(2) doublets ip u = (j'A and ip 2i = ) 

were introduced in Ref. 0. Here we introduce two 
spin-0 boson fields b a , a = 1,2 forming another dou- 
blet bi = (b"J- We then form SU{2) singlets to repre- 
sent the following physical operators fij = \fl li B OL fjfj3i, 
cm = bUu/V2 = {b\j u + b\jl)/V2, and c 2l = 
b\^/V2 = {b\j2i - bljl)/V2. The t-J Hamiltonian 
' Sj - jniUj) - t{& ai c a j + h.c.)] can now be 
written in terms of our fermion-boson (FB) fields. The 
Hilbert space of the FB system is larger than that of 
the t-J model. However, the local SU(2) singlets sat- 
isfying (^ijj^fipai + bj^&i) |phys) = form a subspace 
that is identical to the Hilbert space of the t-J model. 
On a given site, there are only three states that sat- 
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isfy the above constraint. They are /i|0), f 2 \0), and 
+ 62/2 /1) 1 0) corresponding to a spin up and down 
electron, and a vacancy respectively. Furthermore, the 
FB Hamiltonian, as a SU(2) singlet operator, acts within 
the subspace, and has same matrix elements as the t-J 
Hamiltonian. 

Following the standard approach, we obtain the follow- 
ing MF Hamiltonian jllj for the FB system: 

H m = -vJ2 b h - E a o4^^ m + 

i i 

3-7 

E Y ilX ^ 2 + |Aij|2 + toP^oi) + WiUijbj + h.C.) (1) 

(hi) 

where C% = ( ] ). 

The first two terms of H m are included to impose the 
constraints. The value of n is chosen such that the total 
boson density (which is also the density of the holes in 
the t-J model) is (b\bi) — (b^fiu + b' 2i b2i) = x. The val- 
ues of a l 0i are chosen such that {\i> a iT l i>ai + blr l bi) = 0. 
For I — 3 we have 

{fl i f a i + b\ i b u -b\ i b2i) = l (2) 

We see that unlike the U(l) case the density of the 
fermion (/jjj/crf) is not necessarily equal 1 — x. This is 
because a vacancy in the t-J model may be represented 
by an empty site with a b\ boson, or a doubly occupied 
site with a 02 boson. We also notice that the MF Hamil- 
tonian is invariant under local SU (2) transformations, 
W t £ SU(2): -> Wiip ai , h -> Wibi, U l3 WiUyW], 
and a l 0i T l -» W l a\ i T l w\. 

We have searched the minima of the MF free energy for 
the MF ansatz with translation, lattice and spin rotation 
symmetries. We find a phase diagram with six different 
phases (Fig. 1). 

(1) Staggered flux (sF) phase: 

(U i)i+jk = -T 3 x-i(-)^A 
\ViAW = -r 3 X + zR^A {6) 

and a oi = 0. In the U(l) slave-boson theory, the stag- 
gered flux phase breaks translation symmetry. Here 
the breaking of translational invariance is a gauge ar- 
tifact. In fact, a site dependent SU(2) transforma- 
tion Wi — exp(i(— iy*+ l y |n) maps the sF phase to 
the d-wave pairing phase of the fermions: Ui t i+x,y — 
— XT3 ± Ari, which is explicitly translationally invariant. 
In the sF phase the fermion and boson dispersion are 

given by ±Ef and ±E b , where Ef = \J{tj ~ af,) 2 + r)"j, 

e f = (cos k x + cos k y )x, Vf = -^(cos k x - cos k y )A, 
and a similar result for Et with ^ replaced by 2t. Since 
al = we have (flj al ) = 1 and (b\h) = {b\b 2 ) = f . 

(2) The 7r-flux (irF) phase is the same as the sF phase 
except here \ = A. 



(3) The uniform RVB (uRVB) phase is described by 
Eq. (§) with a l 0i = A = 0. 

(4) A localized spin (LS) phase has LZy = and 
a l oi = 0, where the fermions cannot hop. 



U . 3 


"LS / 

- / uRVB 




T/J 




fl; 


. 15 


7 sc 





0.1 0.2 

X t/ J 



FIG. 1. 5(7(2) MF phase diagram for t/J = 1. The 
phase diagram for t/J — 2 is quantitatively very similar to 
the t/J — 1 phase diagram, when plotted in terms of the 
scaled variable xt/J, except the 7rF phase disappears at a 
lower scaled doping concentration. 

(5) The d-wave superconducting (SC) phase is de- 
scribed by U lii+X .y = -XT3 ± Ati and 0% ^ 0, a ' 2 = 0, 
(61) 7^ 0, (62) = 0. Notice that the boson condenses 
in the SC phase despite the fact that in our MF theory 
the interactions between the bosons are ignored. The SC 
MF solution provides an interesting example of finite- 
temperature free boson condensation in two dimensions. 
To see this, notice that the Og term in the FB Hamilto- 
nian - J^i a lUlafai J rb\ i bi l -b\ i b 2 i-l) makes Oq behaves 
like the chemical potential of the fermions. The fermions 
favor a non-zero Og. Let us assume ag < 0, which makes 
(faifai) < 1. A negative Og also makes the 61-band bot- 
tom to be higher than that of b 2 , and the thermally ex- 
cited bosons satisfy {b\bi) therm < (b\b 2 )therm- Thus the 
thermally excited bosons alone cannot satisfy the con- 
straint in Eq.(||). The b\ boson are forced to condense 
at the bottom of the b\ band to satisfy the constraint, 
in the same way that ordinary boson condense to satisfy 
the density constraint. Due to the fermion contribution, 
the total free energy can still be lowered by generating a 
finite ag at low temperatures. 

(6) The Fermi liquid (FL) phase is similar to the SC 
phase except that there is no fermion pairing (A = 0). 

In the following we would like to discuss some sim- 
ple physical properties of the MF phases. uRVB, sF, 7rF 
and LS phases contain no boson condensation and corre- 
spond to unusual metallic states. Since Oq = 0, the area 
of the fermion Fermi sea in the uRVB phase is pinned 
at I of the Brillouin zone. As temperature is lowered, 
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fermions condense in pairs, or develop staggered flux and 
the uRVB phase changes into the sF or irF phases. A 
gap is opened at the Fermi surface which reduces the low 
energy spin excitations. Thus the sF and 7rF phases cor- 
respond to the high temperature spin gap phase. The SC 
phase contains both the boson and the fermion-pair con- 
densations and corresponds to a d-wave superconducting 
state of the electrons. The FL phase containing boson 
condensation corresponds to a Fermi liquid phase of elec- 
trons. However, the area of the Fermi sea produced by 
the SU{2) MF theory is larger than that predicted by 
the Luttinger theorem which reveals a drawback of the 
SU(2) MF theory at high doping concentrations. This is 
probably related to another drawback of the SU(2) the- 
ory that superconducting T c goes down too slowly beyond 
the optimal doping. It appears that the U(l) MF theory 
is better at higher doping. 




FIG. 2. (a) The electron spectral function ImGu in the 
sF phase for t/J = 2,x = 0.041, T/J = 0.13, where % = 0.57 
and A = 0.22. The sharp peaks near uo — are quasi-particle 
peaks. The insert shows a quarter of the Brillouin zone, (b) 
The spectral functions ImGu for three linear scans along the 
line A, B, and C in the insert. 

We would like to point out that the different MF phases 
contain different gauge symmetries. The uRVB and the 
7rF phases have the full SU(2) gauge symmetry. In the 
sF phase the SU (2) gauge symmetry is broken down to 
17(1), (12) while in the SC and FL phases the SU{2) gauge 
symmetry is completely broken. 

Next we calculate the physical electron Green function 
in the sF phase. Using the expression of c a , the MF ap- 
proximation Go is given by the convolution of fermion 
and boson Green functions. The expression of Go is 
lengthy, but can be approximated at low temperatures 
by 



The first term describes the coherent motion of electrons 
with the fermion dispersion. The new feature is the ap- 
pearance of the coherent factors u(k) = *J ~^§^- s g n ( ? 7/) 

and v(k) = \J E ^/ t ' m ^ ne secon d term is the incoher- 
ent background which mainly reflect the boson density 
of states. ImGm exists only for to < and contributes 
1/2 to a total spectral weight (1 + x)/2. 

The coherent part of Go produces only Fermi points 
at S = (±7r/2, ±7r/2). Another feature is that the oc- 
cupied part of the spectral weight of Go contains 1 + f 
electrons as opposed to 1 — x electrons. These unsat- 
isfactory features are due to the absence of correlation 
between fermions and bosons in arriving at Go- In re- 
ality there is a strong attraction between them due to 
gauge fluctuations. In the limit of a single hole, this 
attraction can lead to a bound state with the quantum 
number of a electron, as emphasized in Ref. |lq |. In the 
case of finite hole concentration, we expect that fermion 
particle-hole pairs may be spontaneously excited out of 
the MF ground state so that the 62 (61) bosons can bind 
to fermions (anti-fermions) . The result would be low ly- 
ing physical electron excitations which may resemble a 
Fermi surface. In order to capture this physics, even at a 
very crude level, we assume that after screening, the a l 
fluctuations induce the following short range interaction 
— J2i Ui/j ai T l ipai b\r l bi, where U < 0. We have calcu- 
lated Gu which included this attraction in the Bethe- 
Salpeter approximation. The results are shown in Fig. 
2 and 3. To interpret these results, we consider instead 
G'u = (Go 1 + U)- 1 which we have found numerically 
to be nearly identical to Gu. It is easy to see that G' v 
also has the form G' v — G' coh + G' in . In the case when 
either |u| or \v\ « 1, we find that G' coh = u± i f+llf with 

(If = 2 (\ + UG- ) an< ^ ^ = (l + UG- ) 2 ^ e See a ne S a " 

tive U generates a negative fi / which produces small hole 
pockets. A negative U also enhance the spectral weight 
of unoccupied quasi-particles. This allows us to choose U 
by requiring that the occupied spectral weight in Gu is 
l — x. For general u and v, we find that the coherent part 
of G'u produces pocket-like Fermi surfaces near £ (see Fig 
2) which are determined by 2Ef(l+UGi n ) = Ux(u 2 —v 2 ). 

The quasi-particle weight at the Fermi surface Z = "-^fi 
vanishes at X and is very small on the outer edge of the 
pocket, making it hard to detect. As we approach the 
uRVB phase, A decreases and the pockets are elongated, 
while their area increases with doping. Eventually the 
inner edges of the pockets join together to form a large 
Fermi surface, with low lying excitations near M and a 
shape which resembles the experiment. [0 

The incoherent background of ImGu contains two 
broad peaks separated by 2t, shown in Fig 2a at S and M. 
This follows from the pseudo gap in the boson density of 
states and is a direct consequence of the staggered flux. 
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Coherence factors cause a transfer of spectral weight from 
the low energy to the high energy peak as one goes from 
r to S to X, in qualitative agreement with exact 
diagonalization results. Il5| 




EM r X M 

FIG. 3. The points describe the dispersion of the 
quasi-particle peaks for the sF phase in Fig. 2. The vertical 
bars are proportional to the peak values of ImGu which reflect 
the quasi-particle weight. The solid curve is the quasi-particle 
dispersion for a uRVB phase with the same doping concen- 
tration but at a higher temperature (T/J = 0.19). 



and the full ellipse of the Fermi pocket is completed based 
on our analytic results on G' v . 

The SU(2) MF theory shares many similar physical 
properties with the U(l) MF theory (where the spin gap 
is generated by the d-wave pairing of the fermions) . How- 
ever, there are some qualitative distinctions. 1) The 
d-wave state in the U{\) theory does not produce the 
double-peak structure in the incoherent background of 
ImG. One needs to use a flux phase in the U(l) theory 
to produce the double-peak (at the expense of breaking 
translation or time reversal symmetry Jlif). 2) The d- 
wave state in the U(l) theory does not have Fermi pock- 
ets at finite doping, even if we include the gauge inter- 
action as we did in the SU(2) theory. 3) The sF phase 
in the SU (2) theory contains a gapless 17(1) gauge field, 
which is absent in the corresponding d-wave state in the 
U(l) theory. It has been pointed out that the existence 
of a mass gap in the U(l) theory may de-stablize the d- 
wave state. jl6j The sF phase may be more stable from 
this point of view. 

We would like to thank Z.X. Shen and M. Sigrist for 
very helpful discussions. PAL is supported by NSF- 
MRSEC grant DMR-94-00334 and XGW is supported 
by NSF grant DMR-94-11574 and A.P. Sloan fellowship. 



Fig 3 shows the dispersion of the quasi-particle peaks 
and their spectral weights. Comparing the sF and uRVB 
dispersions, we see clearly a splitting of the spectral 
weight along £ to M which is remarkable for a trans- 
lational invariant state. In our theory the splitting is 
naturally related to the spin gap. We also note that at 
half filling, the bottom of the band at M and T are degen- 
erate. As doping is increased, the spin gap shrinks and 
the occupied band near M moves up in energy, eventually 
producing the flat band near the Fermi surface seen in 
photo-emission experiments. In our calculation, the band 
at r has been pushed up in energy from its bare value 
near — J to — J/2 due to the inclusion of U , whereas the 
experimental value is closer to —2 J. We believe this fea- 
ture, as well as the large spectral weight near T, is an 
artifact of our crude treatment of the gauge fluctuation. 

Fig 2b shows, in more detail, how the Fermi surface 
disappears in the sF phase. In scan C, a ghost band be- 
low the Fermi energy is quite visible after the main peak 
goes above the Fermi surface. This band is connected to 
the occupied band in Fig 3 as the M point is approached. 
The apparence of the ghost band and the sudden reduc- 
tion of the quasiparticle spectral weight in the ghost band 
have been observed in the insulating cuprates. (h]] The 
inner edge of the Fermi pocket in the insert of Fig. 2 is 
determined from the position of the quasiparticle peaks 
at uj = 0. The quasiparticle peaks at the outer edge of 
the Fermi pocket are not visible in our numerical result, 
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